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Abstract 

A novel nonlinear effect of anomalously deep penetration of an external radio 
frequency electric field into a plasma is described. A self-consistent kinetic 
treatment reveals a transition region between the sheath and the plasma. 
Because of the electron velocity modulation in the sheath, bunches in the 
energetic electron density are formed in the transition region adjusted to the 
sheath. The width of the region is of order Vt/uj, where Vr is the electron 
thermal velocity, and to is frequency of the electric field. The presence of 
the electric field in the transition region results in a cooling of the energetic 
electrons and an additional heating of the cold electrons in comparison with 
the case when the transition region is neglected. 

PACS numbers:52. 35. Mw, 52.65Ff, 52.65-y, 52.75-d, 52.80.Pi 

The penetration of the electric field perpendicular to the plasma boundary was studied 
by Landau in the linear approximation [|TJ. He showed that an external electric field with 
amplitude Eq is screened by the plasma electrons in the sheath region in a distance of order 
the Debye length, and reaches a value Eq/e in the plasma, where e is plasma dielectric 
constant. In many practical applications, the value of the external electric field is large: the 
potential drop in the sheath region V s h is typically of order hundreds of Volts and is much 
larger than electron temperature T e , which is of order of a few Volts; and the field penetration 
has to be treated nonlinearly. The asymptotic solution of sheath structure has been studied 
by Lieberman in the limit V s h » T e [Q. In this treatment, the plasma sheath boundary 
is considered to be infinitely thin and the position of the boundary is determined by the 
condition that the external electric field is screened in the sheath regions when electrons 
are absent. Electron interactions with the sheath electric field are traditionally treated as 
collisions with a moving potential barrier (wall). It is well known that multiple electron 
collisions with an oscillating wall result in electron heating, provided there is sufficient 
phase-space randomization in the plasma bulk. It is common to describe the sheath heating 
by considering the electrons as test particles, and neglecting the plasma electric fields J|. 
Kaganovich and Tsendin proved in Ref. j|] that accounting for the electric field in the plasma 
reduces the electron sheath heating, and the electron sheath heating vanishes completely in 
the limit of uniform plasma density. Therefore, an accurate description of the rf fields in the 
bulk of the plasma is necessary for calculating the sheath heating. The electron velocity is 
oscillatory in the sheath, and as a result of this velocity modulation electron density bunches 
appear in the region adjusted to the sheath. The electron density perturbations decay due 
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to phase mixing over a length of order Vt/uj, where Vr is the electron thermal velocity, 
and uj is the frequency of the electric field. The electron density perturbations polarize the 
plasma and produce an electric field in the plasma bulk. This electric field, in turn, changes 
the velocity modulation and correspondingly influences the electron density perturbations. 
Therefore, electron sheath heating has to be studied in a self-consistent nonlocal manner 
assuming a finite temperature plasma. 

Notwithstanding the fact, that particle- in-cell simulations results are widely available for 
the past decade [5-7] a basic understanding of the electron sheath heating is incomplete, 
because no one has studied the electric field in the plasma bulk using a nonlocal approach, 
similar to the anomalous skin effect for inductive electric field ||. In this regard, analytical 
models are of great importance because they shed light on the most complicated features 
of collisionless electron interactions with the sheath. In this Letter, an analytical model is 
developed to explore the effects associated with the self-consistent non-local nature of the 
phenomenon. 

One of the approaches to study electron sheath heating is based on a fluid description of 
the electron dynamics. For the collisionless case, closure assumptions for the viscosity and 
heat fluxes are necessary. In most cases, the closure assumptions are made empirically or 
phenomenologically [6, 7]. The closure assumptions have to be justified by direct comparison 
with the results of kinetic calculations as is done, for example, in Ref. 0. Otherwise, 
inaccurate closure assumptions may lead to misleading results as discussed below. 

To model the sheath-plasma interaction analytically, the following simplifying assump- 
tions have been adopted. The discharge frequency is assumed to be small compared with 
the electron plasma frequency. Therefore, most of the external electric field is screened in 
the sheath region by an ion space charge. The ion response time is typically larger than the 
inverse discharge frequency, and the ion density profile is quasi-stationary. There is an ion 
flow from plasma bulk towards electrodes. In the sheath region, ions are being accelerated 
towards the electrode by the large sheath electric field, and, the ion density in the sheath 
region is small compared with the bulk ion density. In the present treatment, the ion density 
profile is assumed fixed and is modeled in a two-step approximation: the ion density is 
uniform in the plasma bulk, and the ion density in the sheath n s h < is also uniform 
(see Fig.l). At the sheath-plasma boundary, there is a stationary potential barrier for the 
electrons (e§ s h), so that only the energetic electrons reach the sheath region. The potential 
barrier is determined by the quasineutrality condition, i.e., when the energetic electrons 
enter the sheath region, their mean density is equal to the ion density [n e (§ s h) = n sh \. 

The electron density profile is time-dependent in response to the time-varying sheath 
electric field. The large sheath electric field does not penetrate into the plasma bulk. There- 
fore, the quasineutrality condition holds in the plasma bulk, i.e., the electron density is equal 
to ion density, n e = n\>. In the sheath region, the electrons are reflected by the large sheath 
electric field. Therefore, n e = n s h for x > x s h(t), and n e = for x < x s h(t), where x s h(t) is 
the position of the plasma-sheath boundary ||. From Maxwell's equations it follows that 
V ■ J = 0, where the total current J is the sum of the displacement current and the electron 
current. In the one-dimensional case, the condition V • J = yields the conservation of the 
total current: 
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where jo is the amplitude of the rf current controlled by the external circuit and is the 
initial phase. In the sheath, electrons are absent in the region of large electric field, and the 
Eq.([L]) can be integrated to give || 

E sh (x,t) = — — [— 1 - cos(u;t + </>)] + 4ir\e\n sh x, x < x sh (t) (2) 

where Poisson's equation has been used to determined the spatial dependence of the sheath 
electric field. The first term on the right-hand side of Eq.(^j) describes the electric field at 
the electrode, the second term relates to ion space charge screening of the sheath electric 
field. The position of the plasma-sheath boundary x s h{t) is determined by the zero of the 
sheath electric field, E sh [x s h(t),t] = 0. From Eq.(^) it follows that 

x afc (t) = — [l + cos(urt + </>)], (3) 

where V s ho = jo/(en sh ) is the amplitude of the plasma-sheath boundary velocity. The ion 
flux on the electrode is small compared with the electron thermal flux. Because electrons 
attach to the electrode, the electrode surface charges negatively, so that in a steady-state 
discharge, the electric field at the electrode is always negative, preventing an electron flux on 
the electrode. However, for a very short time (u)t n + <p ~ 7r(l + 2n)) the sheath electric field 
vanishes, allowing electrons to flow to the electrode for compensation of the ion flux. Note 
that there is a large difference between the sheath structure in the discharge and the sheath 
for obliquely incident waves interacting with a plasma slab without any bounding walls. 
Because electrodes are absent, electrons can move outside the plasma, and the electric field 
in the vacuum region, E s h(x, t) = (4Trj Q /u) cos(ut + <f>), may have a different sign. Therefore, 
electrons may penetrate into the region of large electric field during time when E s h(x, t) < 
[10,11]. However, in the discharge, because the sheath electric field given by Eq.fl2|) is 
always reflecting electrons, the electrons never enter the region of the large sheath electric 
field, which is opposite to the case of obliquely incident waves. 

The calculations based on the two-step ion density profile model is known to yield dis- 
charge characteristics in good agreement with experimental data and full-scale simulations 
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Throughout this paper, linear theory is used because the plasma-sheath boundary ve- 
locity and the mean electron flow velocity are small compared with the electron thermal 
velocity [4,5]. The important spatial scale is the length scale for phase mixing, l mix = Vt/oj. 
The sheath width satisfies 2V s ho/uj « l mix because V s n « Vt- Therefore, the sheath width 
is neglected, and electron interactions with the sheath electric field are treated as a bound- 
ary condition. The collision frequency (u) is assumed to be much less than the discharge 
frequency [y << u), and correspondingly the mean free path is much larger than the length 
scale for phase mixing. Therefore, the electron dynamics is assumed to be collisionless. The 
discharge gap is considered to be sufficiently large compared with the electron mean free 
path, so that the influence of the opposite sheath is neglected. The effects of finite gap width 



are discussed in Ref. [13 



The electron interaction with the large electric field in the sheath is modelled as collisions 
with a moving oscillating rigid barrier with velocity V s h{t) = dx s h(t)/dt. An electron with 
initial velocity — u after a collision with the plasma-sheath boundary - modeled as a rigid 
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barrier moving with velocity V s h{t) - acquires a velocity u + 2V s h- Therefore, the power 
deposition density transfer from the oscillating plasma-sheath boundary is given by [§] 



m 

sh = l> 



IT du[u + V sh (t)]\(2V sh (t)+u) 2 -u 2 } f sh (-u,t)\, (4) 



where m is the electron mass, f s h(—u,t) is the electron velocity distribution function in 
the sheath, and (• • •) denotes a time average over the discharge period. Introducing a new 
velocity distribution function g(—u',t) = f s h[—u — V s h(t),t], Eq.(f|) yields 

Psh = ~2m (v sh (t) J™ u' 2 g(-u', t)du'\ , (5) 

where —u' = —u — V s h is the electron velocity relative to the oscillating rigid barrier. From 
Eq.(^|) it follows that, if the function g(u') is stationary, then {P s h = 0) there is no collisionless 
power deposition due to electron interaction with the sheath [7, 14]. For example, in the limit 
of a uniform ion density profile n s h = n&, qW) is stationary {in an oscillating reference frame 
of the plasma- sheath boundary), and the electron heating vanishes H|. Indeed, in the plasma 
bulk the displacement current is small compared with the electron current, and from Eq.([I|) it 
follows that the electron mean flow velocity in the plasma bulk, Vb(t) = —jo sin(o;i+0)/|e|n&, 
is equal to the plasma-sheath velocity V s h(t), from Eq.(D). Therefore, the electron motion 
in the plasma is strongly correlated with the plasma-sheath boundary motion. From the 
electron momentum equation it follows that there is an electric field, E b = m/edVb(t)/dt, 
in the plasma bulk. In a frame of reference moving with the electron mean flow velocity, 
the sheath barrier is stationary, and there is no force acting on the electrons, because the 
electric field is compensated by the inertial force (eEb — mdVb(t)/dt = 0). Therefore, electron 
interaction with the sheath electric field is totally compensated by the influence of the bulk 
electric field, and the collisionless heating vanishes ||. 

The example of a uniform density profile shows the importance of a self-consistent treat- 
ment of the collisionless heating in the plasma. If the function g(u',t) is nonstationary, 
there is net power deposition. In this Letter, a kinetic calculation is performed to yield the 
correct electron velocity distribution function g(u',t) and, correspondingly, the net power 
deposition. 

The electron motion is different for the low energy electrons with initial velocity in the 
plasma bulk \u\ < u s h, where u 2 sh = 2e& s h/m, and for the energetic electrons with velocity 
\u\ > u s h. The low energy electrons with initial velocity in the plasma bulk — u are reflected 
from the stationary potential barrier e$ s ^, and then return to the plasma bulk with velocity 
u. High energy electrons enter the sheath region with velocity u\ = —(u 2 — u 2 h Y^ 2 . They 
have velocity u<i = 2V s h — U\ colliding with the moving rigid barrier, and then return to the 
plasma bulk with velocity (u 2 , + u 2 sh ) l l 2 |TJ . 



As the electron velocity is modulated in time during reflections from the plasma-sheath 
boundary, so is the energetic electron density (by continuity of electron flux). This phe- 



nomenon is identical to the mechanism for klystron operation fl6| . The perturbations in 
the energetic electron density yield an electric field in the transition region adjusted to the 
sheath. 

The electron velocity distribution function is taken to be a sum of a stationary isotropic 
part fo(u) and a nonstationary anisotropic part fi(x, u, t). /i is to be of the form f±(x, u, t) = 
fi(x, u) exp(-iu)t) . The linearized Vlasov equation becomes 
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(6) 



where the term on the right-hand side accounts for rare collisions {y << uj). All 
time-dependent variables are assumed to be harmonic functions of time, proportional to 
exp(— iut), and, in the subsequent analysis, the multiplicative factor exp(— iut) is omitted 
from the equations. The electron velocity distribution function must satisfy the bound- 
ary condition at the plasma-sheath boundary (x = 0) corresponding to f(0,u) = /(0, —u) 
for \u\ < u sh , and f s h{u') = f s h(2V sh - u'), for u > u sh , where v! = (u 2 - u 2 sh ) 1/2 and 
f s h is the electron velocity distribution in the sheath. From energy and flux conservation, 
v! f s h(u')du' = uf(u)du, it follows that f s h(u') = f[(u' 2 + u 2 sh ) l l 2 }. Linearly approximating 
the boundary conditions yields 

/i(0,u) = /i(0,-u), 0<u<u sh , (7) 



/ 1 (0,iO = A(0,-iO+2V, h -f^, u>u sh . (8) 

u du 

The electric field is determined from the condition of conservation of the total current (jo), 
which gives 

/OO jjjj 
ufi(x,u)du - —E(x) = j, (9) 
-oo 47T 

where j = jo^ l ^ +n ^ 2 \ and the first term is the electron current and the second term corre- 
sponds to a small displacement current. Equations (Q) and @, together with the boundary 
conditions (0), (fD comprise the full system of equations for the bulk plasma. 

It is convenient to solve Eq. (|6|) by continuation into the region x < 0. First, we introduce 
the artificial force 

F(x, u) = 2mV sh u'5(x)e(\u\ - u sh ), (10) 

where V s h = j/en s h, S(x) is the Dirac delta-function, and Q(u) is the Heaviside step function. 
The force in Eq. fllQ|) accounts for the change of the energetic electron velocity in the sheath 
region. Equation (|6|) together with the boundary conditions (0) and (§) are equivalent to 
Eq. (|jD with the force in Eq. fllCf ) added to the third term of Eq. @ . This gives 

... dfx eE(x) + F(x,u)df 

-iujf 1 + u— + — = -vf 1} (11) 

ox m du 

where the boundary condition ([/]) for all electrons becomes 

/ 1 (0,tt)=/ 1 (0,-«). (12) 

In this formulation, the half-space problem is equivalent to that of an infinite medium 
in which the electric field is antisymmetric about the plane x — 0, with E(x) = —E(—x) 
[1, 17]. Such a continuation makes Eq. ( |TTD invariant with respect to the transformation 
x — > —x, and u — > —u. Electrons reflected from the boundary in the half-space (x > 0) 
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problem correspond to electrons passing freely through the plane x = from the side x < 
in the infinite-medium problem. 

A spatial Fourier the transform of Eq. ( |TTD gives 

eE(k) + F sh (u) df 

= —7 , , . s ~r, (13) 

rm(u; — uk + z^J aw 
where is the Fourier transform of E(x) 

/oo 
E(x)exp(-ikx)dx, (14) 
-oo 

and F s h(u) = 2mV s hu'Q(\u\ — u s h). It is convenient to divide the electric field in the plasma 
into two parts corresponding to E(x) = Ei(x) + E b sgn(x), where Ei(x) — > for x — > oo, and 
E b is the value of the electric field far away from the sheath region. The Fourier transform 
of the electron current can be obtained by integrating Eq. ( |T3|) over velocity, yielding 

j(k) = a{k)E 1 {k) - j [E sh a sh (k) + a(k)E b ] , (15) 



ie 2 Z" 00 u df 



a(k) = / t : ——du, (16) 

m J-oo {oo — uk + iv) du 

ie 2 k roo uu'Q(\u\ - u sh ) df 

(Jsh(k) = r — — / . — du, (17) 

(uj + %v)m J-oo (uj — uk + iv) du 

where a(k) is the electron conductivity, cr s h(k) is the effective conductivity due to elec- 
tron interaction with the sheath, and E s h = {—iuo + v)mV s f h /e is the effective electric field 
corresponding to V s h- 

The Fourier amplitude E\{k) is to be determined from Eq.(^|) continued into the half- 
space x < 0. Because E(x) is an antisymmetric function about the plane x = 0, jo is 
continued with negative sign into the half-space x < 0, and the Fourier transform of josgn(x) 
is —2ij /k. Substituting E{k) = E 1 (k) — 2iE b /k and j = [cr(0) — iu/Air] E b into Fourier 
transform of Eq. (||) gives 

E (k) = 2^ [a(0)-a(k)}E b -E sh a sh (k) 
U ' k a(k) - g 1 ' 

Notice that, if the plasma density in the sheath is equal to the bulk density n s h = n b , then 
u s h — 0, E b — E s h and cr(0) — a{k) = a s h{k). Therefore, E\{k) = and the uniform electric 
field E b satisfies the current conservation condition, as discussed earlier. 
The profile for E\(x) given by inverse Fourier transform 

1 r°° 

E(x) = — E{k)exp{ikx)dk (19) 
2n J-oo 

is shown at the top in Fig.2. For x < 6Vr/ uj the electric field profile is close to E\(x) ~ 
£i(0) exp(-\xio/V T ), where Ex(0) = -0.72, and A = 0.19 + 0.77i for the conditions in 
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Fig. 2. For x > 6Vr/ uj, the electric field profile is no longer a simple exponential function, 
similar to the case of the anomalous skin effect [|l7j]. The three components of current 
corresponding to the first, second, and third terms in Eq. flT5| ) are shown at the bottom 
in Fig. 2. The first term describes the current (j tr ) driven by the electric field Ex(x) under 
the assumption of specular reflection at the boundary. The second term relates the current 
(jsh) of the energetic electrons owing only to a velocity change due to reflections from 
the large sheath electric field. The third term describes the current (jb) driven by the 
uniform electric field Eb under the assumption of specular reflection at the boundary. Due 
to the boundary condition of specular reflection in Eq. (0), both of the currents jb and 
jtr are equal to zero at x — 0. Also, both of the currents j tr and j s h vanish at x > 
15Vt/ uj due to phase mixing, and the only current left here is jb- In contrast to large x, 
at small x « Vt/ uj the total current is entirely due to energetic electrons interacting 
with the sheath j s h. Indeed, the energetic electrons enter the sheath region with velocity 
distribution f s h{u'). The electron current is given by the sum of the contribution from the 
electrons approaching the oscillating barrier and from the electrons already reflected from 
the barrier, j sh = J^ h u' f sh (u')du' + jZ^ h u' f sh (u')du' . Because f sh {u') = f s h(2V sh - u'), 
j sh = 2eV sh jy sh f sh (u)du « eV sh f sh (u')du' = en sh V sh = j a sm(ut + 4>). In the last 
calculation the contribution to the density by electrons with velocity u < V^o is omitted. 
Their contributions are second-order effects in V s ho/Vr, which are neglected in the present 



study |[L5| . Therefore, in the sheath region, when electrons are present, and in the nearest 
vicinity of the sheath all current is conducted by the energetic electrons. As can be seen 
in Fig.2, the current conservation condition, jt r (x) + jsh(x) + jb(x) = [c(0) — iuj/Air} Eb, is 
satisfied for arbitrary x. 

The difference in phase of the currents of the energetic and low energy electrons was 
observed in Ref. ||, but it was misinterpreted as the generation of electron acoustic waves. 
Electron acoustic waves can be excited if the denominator of the right-hand side of Eq. ([18]) 
has a pole at frequency uj, which corresponds to the root of the plasma dielectric function, 
e = 1 + Airier (k)/uj. For a Maxwellian electron distribution function, the pole does not exist 



for uj « uj p , where u p = y4ire 2 rib/m is the electron plasma frequency. But the electron 
acoustic waves can exist if the plasma contains two groups of electrons having very different 
temperatures The wave phase velocity is u/k = \Jn c /nh\jTh/m , where n c and rih are 
the electron density of cold and hot electrons, respectively, and Th is the temperature of the 
hot electrons. The electron acoustic waves are strongly damped by the hot electrons, unless 
n c « n h and T c «T h , where T c is the electron temperature of the cold electrons In 
the opposite limit, n c > 4n h , the electron acoustic waves do not exist |Tj|. In capacitively- 
coupled discharges, the electron population does stratify into two populations of cold and 
hot electrons, as has been observed in experiments and simulation studies [19,20]. Cold 
electrons trapped in the discharge center by the plasma potential do not interact with the 
large electric fields in the sheath region and have low temperature. Moreover, because of 
the nonlinear evolution of plasma profiles, the cold electron density is much larger than the 
hot electron density p0| . Therefore, weakly-damped electron acoustic waves do not exist 



in the plasma of capacitively-coupled discharges. Reference used the fluid equation and 
neglected the effect of collisionless dissipation, thus arriving at the wrong conclusion about 
the existence of weakly-damped electron acoustic waves. 

The power deposition is given by the sum of the power transferred to the electrons by 
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the oscillating rigid barrier in the sheath region and by the electric field in the transition 
region, 

Ptot = Psh + Ptr. (20) 

Here P s h is given by Eq.(|]), which after linearization yields 

Psh = PshO + Pshl. (21) 



In Eq. (|21|) , P s ho is the power dissipation in the sheath neglecting any influence of electric 
field, 



shO 



2m 2u' V sh {tf foshi-u', t)\ , (22) 



and P s hi accounts for the influence of the electric field on fi and correspondingly on the 
power dissipation in the sheath, 



shl 



2m / jH V sh (t)u' 2 hsh{-u\ x = 0, t)du'\ . (23) 



Time averaging, changing variables from v! to u, and integration by parts in the first term 
yield 

P sh = mj™ \-\V sh \ 2 u' 2 ^ + Re [V sh u'uf*(-u, x = 0)]| Q(\u\ - u sh )du, (24) 
where /j* is solution to Eq.(H), 

f*(-u, x = 0) = — [°° E\x)e-^ +u)x/u dx. (25) 



mu du Jo 

Time averaging the power deposition in the transition region, J °° (jE) dx, gives 

1 f°° 

P tr = - Re / j E*dx. (26) 

2 jo 

Substituting j = ^ n bVb, where = eE b /muj is the amplitude of the mean electron flow 
velocity in the plasma bulk and = was assumed in Eq.(Q), we obtain P^ = 1/2 KejoE* = 
— l/2eribVb Im E\(x). Therefore, Pb is determined by the imaginary part of E%, and can 
be either positive or negative (see Fig. 2). Negative power density has been observed in 
numerical simulations |J. 

Substituting j = ]e + j s h, where ]e = jb + jtr, the power deposited by the current je 



can be calculated by continuing into infinite space and using the Fourier transform |L7 



-Re/ j E E*dx = -Re j E E*dx = —Re j E (k)E*(k)dk, (27) 

2 JO 4 J-oo 07T J-oo 

where jE{k) = cr(k)E(k). Finally, substituting the conductivity from Eq.(16) yields 

iRer^=-ir*zii^„. (28) 

2 Jo 4 jo m du 



The current j s h is determined by the perturbed electron velocity distribution function due to 
reflections from the sheath electric field. The perturbed distribution function / ls /j at x = 
is given by Eq.(||), and for x > the solution to the Vlasov equation becomes 



f lsh (x,u) = -2V sh -^e("W\ 
u du 



(29) 



Calculating the current j s h by integrating fuh from Eq.(p9"D over velocity, and substituting 
the current into Eq.(|26|) gives 



j />0O /■OO 

-Re / j ah E*dx = - Re V sh / u'E*(k 
2 Jo Jo 



u du 



du 



(30) 



Substituting /* from Eq. Q2"5D into Eq. 
and (|30|) yield 



Ptot = ~ [ muD u (u)—^-du 
Jo du 

where D u {u) is the diffusion coefficient in velocity space, 

u | du | 2 



and adding the contributions from Eqs. 



(31) 



DM 
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(32) 



and du is the change in electron velocity after passing through the transition and sheath 
regions, 



du = 2iVh 



u rib 
u n sh 



6(|w| - u sh ) - 1 



eEi(k = to/u) 



u 



(33) 



A plot of \du\ 2 /2 is shown in Fig. 3. Taking into account the electric field in the plasma (both 
E b and E\) reduces \du\ for energetic electrons (u > u s h) and increase \du\ for slow electrons 
(u < u s h). Therefore, the electric field in the the plasma cools the energetic electrons and 
heats the low energy electrons, respectively. Similar observations were made in numerical 
simulations 0. 

Figure 4 shows the dimensionless power density as a function of n^/rtsh. Taking into 
account the electric field in the plasma (both E b and Ei) reduces the total power deposited 
in the sheath region. Interestingly, taking into account only the uniform electric field Ej, 
gives a result close to the case when both Eb and E\ are accounted for. The electric field E\ 
redistributes the power deposition from the energetic electrons to the low energy electrons, 
but does not change the total power deposition (compare Fig.3 and Fig. 4). Therefore, the 
total power deposition due to sheath heating can be calculated approximately from Eq. ( PH ) . 
taking into account only the electric field E b . This gives 



Pi 



tot 



-mVS 



u 



u rib 
u n sh 



-\ 2 



6(it - u sh ) 



dfo 
du 



du. 



(34) 



The result of the self-consistent calculation of the power dissipation in Eq.(|34]) differs from 
the non-self-consistent estimate in Eq . (|22|) by the last term in Eq.([34j), which contributes 
corrections of order n s h/rib to the main term. 
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APPENDIX: 



I. PROPERTIES OF E X {K) 

The Fourier transform E\(k) has the following properties in the limits of small and large 
k. At small k (k << u/Vr), Ei(k) ~ k, because the numerator in the last factor on the 
right-hand side of Eq.QT|)~ k 2 ( [cx(0) - a(k)\ ~ k 2 and a sh (k) ~ k 2 ). Because E x (k) ~ k 



for small k, J Ei(x)dx = similarly to the case of anomalous skin effect ]T7 . 

At large k (r^ 1 >> k » lu/Vt), Ei(k)~l/k, because both the numerator and the 
denominator in the last factor on the right-hand side of Eq.(|l8D are reciprocal to k~ 2 
(a(0)Eb = E s h<J s h(k — > oo)). E\{x) at small x is determined by behavior of E\{k) at 
large k. In the limit of large k (r^ 1 >> k >> uj/Vt) 



where 



Ei(k) = ^, (Al) 



A = E b - ^-E sh . (A2) 



Here, 

B = lim a(k)k 2 , (A3) 
C = lim [a sh (oo) - a sh (k)} k 2 . (A4) 

k— too 

For a Maxwellian electron distribution function, substituting definitions of conductivities 
Eqs.flTEl) and (|I7[) into Eqs.(A3) and (A4), respectively, yields 

B= i£»r» iJA du = ( A5) 

m j-oo u du T 

-ie 2 u [°° u' df sh , -ie 2 um sh , r „ 2 x 

C = / —5- — T~T~f du = m 1 ~ V™ nsh e n.h 1 - erf (u nsh )\ , (A6) 

m J -00 u' 2 + u sh du' T v 7 

where u ns h = u s h/Vr, and eii{u ns h) is the error function. Form Eq.(|l9]), E\(x) at small x is 
given by 

m / r\\ f°° sin(fcx) „ . . . . 

E 1 ( x ^0) = / — ) — -dk = -A. (A6) 

n Jo k 

Substituting and E s h = E^n^jrish and values of B and C from Eqs. (A5) and (A6) into 
Eq.(A2) gives 

£i(0) = -Vnu nsh e u ™h [1 - eri(u nsh )} 
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FIG. 1. Schematic of a sheath. The negatively charged electrode pushes electrons away by 

different distances depending on the strength of the electric field at the electrode. Shown are the 

density and potential profiles at two different times. The solid line is at the time of maximum 

sheath expansion. 
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FIG. 2. Plots of the electric field and the current normalized to their respective values in the 
plasma bulk, and e 2 nEi ) /muJ, as functions of the normalized coordinate xoj/Vt for the following 
parameters: n^/n;, = 1/3, uj/lo p = 1/100, and a Maxwellian electron distribution function. The 
upper graph shows profiles of E\[x): (a) amplitude - solid line; (b) real part - dashed line; (c) 
imaginary part - dotted line; and (d) phase with respect to phase of E\, divided by n- dash-dotted 
line. The lower graph shows profiles of imaginary part of currents: (e) jtr - solid line; (f) jsh 
-dashed line; and (g) jb - dotted line. 




FIG. 3. Plot of the average square of the dimensionless velocity kick as a function of the 
dimensionless velocity for the conditions in Fig.l, taking into account (a) both E\(x) and E\>- solid 
line; (b) only E\>- dashed line; and (c) no electric field - dotted line. 
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FIG. 4. Plot of the dimensionless power density as a function of the ratio of the bulk plasma 
density to the sheath density, taking into account (a) both E\{x) and Ef,- solid line; (b) only Ef,- 
dashed line; and (c) no electric field - dotted line. 
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